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Abstract 



We report a new measurement of J/ip, if}' and Drell-Yan cross-sections, in the 
kinematical domain —0.425 < y cm < 0.575 and —0.5 < coshes < 0.5, performed 
at the CERN-SPS using 400 GeV/c incident protons on Be, Al, Cu, Ag, W and 
Pb targets. The dependence of the charmonia production cross-sections on the size 
of the target nucleus allows to quantify the so-called normal nuclear absorption. 
In the framework of the Glauber model, this new measurement is combined with 
results previously obtained with the same apparatus, under different experimental 
conditions, and leads to a precise determination of the J /if} and if>' absorption cross- 
sections in the surrounding nuclear matter. 



Accepted for publication in Euro. Phys. J. C. 



Universita del Piemonte Orientale, Alessandria and INFN- Torino, Italy 

2 LPC, Univ. Blaise Pascal and CNRS-IN2P3, Aubiere, France 

3 IFA, Bucharest, Romania 

4 Universita di Cagliari/INFN, Cagliari, Italy 

5 CERN, Geneva, Switzerland 

6 LIP, Lisbon, Portugal 

7 INR, Moscow, Russia 

8 IPN, Univ. de Paris-Sud and CNRS-IN2P3, Orsay, France 

9 LLR, Ecole Polytechnique and CNRS-IN2P3, Palaiseau, France 

10 Universita di Torino/INFN, Torino, Italy 

11 IPN, Univ. Claude Bernard Lyon-I and CNRS-IN2P3, Villeurbanne, France 

12 YerPhI, Yerevan, Armenia 

a) also at 1ST, Universidade Tecnica de Lisboa, Lisbon, Portugal 

b) also at Faculty of Physics and Nuclear Techniques, AGH University of Science and 
Technology, Cracow, Poland 

c) on leave of absence of YerPhI, Yerevan, Armenia 

d) deceased 



2 



1 Introduction 



Charmonium production in proton-nucleus collisions has already been extensively 
studied. Several fixed target experiments (NA3 [lj, E772 0, NA38 [3], E866 @], 
E672/E706 [5] and NA50 [6j [7j) have studied and ip' production, including the 
nuclear dependence of their production cross-sections, in the energy range yfs ~ 20- 
40 GeV and as a function of several kinematical variables: Feynman x (xf), rapidity 
(y) and transverse momentum (pt). These measurements have been analysed in view 
of understanding the physics mechanisms underlying charmonia production [8j [9] . 

The study of J/ip production in proton-nucleus collisions is also crucial for a 
correct interpretation of the J/ip suppression patterns experimentally observed in 
heavy-ion collisions, at the CERN/SPS [TO], predicted 20 years ago to signal Quark 
Gluon Plasma formation [11]. It is, indeed, very important to establish a robust base- 
line on the basis of proton-nucleus data, the so-called "normal nuclear absorption", 
with respect to which we can clearly identify new phenomena, specific of high-energy 
heavy-ion collisions. 

This work reports new measurements on J/ip and ip' production from proton- 
nucleus data collected by the NA50 experiment with 400 GeV protons incident on 
various targets. Together with previous proton-nucleus results from the same ex- 
periment at 450 GeV, this new determination leads to a good characterization of 
the charmonia survival probability in "cold nuclear matter" and precisely reveals the 
differences between the two charmonium states, and ip'- 

2 Experimental apparatus and data selection 

The data samples analysed in this paper have been collected with the muon spec- 
trometer used by the NA50 experiment since 1994, as described in [TO]. The target 
region of the experiment, however, has been adapted to the use of six different tar- 
get materials, sequentially exposed to the beam in relatively short time intervals, as 
shortly summarized hereafter. 

The six different targets chosen for these measurements were made of Be, Al, 
Cu, Ag, W and Pb, so as to cover a wide range of nuclear sizes. They were disks 
of 12 mm diameter, much larger than the proton beam spot of 500 /im r.m.s., en- 
suring that no incoming proton misses the target transverse area. Their thicknesses 
(L tgt ) were chosen in order to yield event samples of similar statistics for each target 
species, in data collection periods ( "runs" ) of around 1-2 hours each, while avoiding 
too thick targets (more than 13 cm or 40% of an interaction length), which could 
degrade the measurements. Table [1] summarizes some target properties, the num- 
bers in parentheses giving the relative errors associated with each quantity. The 
target densities, p, were measured, while the interaction lengths, Ai nt , were computed 
using the cross-sections published in [12]. The target effective length is defined as 
L cS = A int [1 - exp (-pL tgt /A int )]. 

A remotely operated rotating target holder placed each different target in the 
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Table 1: Properties of the targets used to collect the data samples analysed in the 
present work. 



beam sequentially, with a target change every hour or so. This feature allowed us 
to balance time-dependent effects among different exposures, reducing to a negligible 
level the systematic uncertainties related to beam characteristics. 

The hadron absorber separating the target region from the muon spectrometer 
was made of 60 cm of Al 2 03, followed by 400 cm of carbon and 80 cm of iron. 
The overall beam time allocated to the measurements reported in this paper was 
only 4 days, forcing us to run with a high intensity proton beam (3-4 x 10 9 protons 
per 3.2 s burst). Nevertheless, the data acquisition worked with average lifetimes 
between 99 % and 96 %, depending on the specific trigger rate of each target exposure. 
The time stability of the targeting efficiency was continuously monitored by three 
scintillator telescopes pointing to the target, with axes orientation at 90° relative 
to the beam line, one vertical and two horizontal. The intensity of the incident 
proton beam was redundantly measured, for continuous cross-checking purposes, by 
three independent ionization chambers placed upstream of the target, with a sensitive 
transverse size which ensured that they intercepted all the incoming beam. The 
ionization current is integrated during every burst and the total collected charge is 
proportional to the number of crossing protons. These detectors have been precisely 
calibrated in a dedicated low intensity data collection using a coincidence of two 
scintillators equipped with very fast electronics. Their linearity was checked up to 
the highest beam intensity used in the experiment and a 3% systematic uncertainty 
on the number of incident protons has been estimated [T3] . 

Two measurements of the dimuon trigger efficiency, made in dedicated runs at the 
beginning and middle of the data taking period, gave (86.3±0.3) % and (88.0±0.3) %, 
respectively. From the study of several different observables as a function of time, a 
2.3% systematic error has been deduced. 

The reconstruction efficiency depends on the occupancy of the spectrometer Multi 
Wire Proportional Chambers which is a function of the luminosity of each exposure 
(target length and beam intensity). It was determined, with its corresponding error, 
from the individual efficiencies of each of the two sets of 12 wire planes detecting 
uncurved tracks, upstream and downstream from the magnet, and cross-checked with 
a Monte-Carlo study in which generated dimuons embedded into real events are 
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processed through the same reconstruction algorithm as applied to real data [H]. 
The final absolute experimental values vary between 77% and 84%, with a 0.2% 
systematic uncertainty, depending on the individual short time exposure. 

The dimuon events are selected according to the standard quality criterium used 
in NA50, mostly driven by geometrical arguments. In particular, each single muon 
must cross the air core toroidal magnet (ACM) of the spectrometer in one of the 
six air-sectors, and it must point to the target. The latter condition is assured by a 
cut on the p x -D t ar g et variable, where p is the momentum of the muon and -Dtarget 
is the transverse distance, in the plane containing the centre of the target, between 
the beam line and the intercept of the extrapolated muon track. Phase space regions 
where the acceptance is less than 10 % of its maximum value are removed by applying 
the following kinematical cuts: —0.425 < y cm < 0.575 and —0.5 < coshes < 0.5, 
where y cm is the dimuon center of mass rapidity and #cs is the polar angle in the 
Collins-Soper frame. 

A small fraction of dimuons from collisions upstream or downstream of the tar- 
get survives the selection criteria. They are estimated from data samples collected 
with no target in place, after proper normalisation to the luminosity integrated in 
the standard runs. The contribution of such spurious dimuons in the mass region 
2.7 < M < 3.5 GeV/c 2 is less than ~ 3 % for the lightest target and less than ~ 1 % 
for the heaviest one. In spite of their small number, they must be carefully taken 
into account because the reconstruction program assumes that these dimuons are 
produced in the target and, therefore, their reconstructed invariant mass is shifted to 
regions where they may contribute more significantly to the total number of events 
(in particular at around 2.8 GeV/c 2 , as seen in Fig. [1] below). 
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Table 2: Some features of the collected data samples. 

Table [2] presents, for each data sample, the total number of incident protons, the 
average efficiency (affected by a 2.3% systematic uncertainty), and the number of 
reconstructed opposite-sign muon pairs in the mass region 2.7 < M < 3.5 GeV/c 2 
after the selection cuts. 
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3 Data analysis 

The opposite-sign muon pair invariant mass distributions are analysed as a superpo- 
sition of several contributions. For masses lower than 2 GeV/c 2 , the mass spectrum 
is dominated by uncorrelated decays of pions and kaons. This background can be 
estimated from the number of measured like-sign pairs, according to the standard 
formula: 



where ACM+ and ACM— denote the magnetic field polarities used to collect the 
events. This relation is valid if the probability to detect a muon in the apparatus (ac- 
ceptance) is independent of its charge. An additional offline cut is therefore imposed 
which rejects events where one (or both) of the muon tracks would not be accepted 
if it had the opposite charge (or if the magnetic field had the opposite polarity). i?b g 
accounts for charge correlations. It tends to unity in events of very high charged 
track multiplicities, where such correlations can be neglected. In the case of our data 
samples, i?b g varies between 1.03 and 1.05, depending on the target, as estimated 
by fitting the mass region 1.5-2.3 GeV/c 2 , where the combinatorial background is 
dominant. These values are very close to unity because the beam intensities we have 
used are so high that it frequently happens that the two muons of a given triggered 
event come from a different p-A collision and are, therefore, uncorrelated. 

For dimuon masses above 2.7 GeV/c 2 , the combinatorial background contributes 
less than 1 % of the total number of events. In this mass region the relevant contri- 
butions are Drell-Yan dimuons and the J /if} and if)' resonances. We have also taken 
into account the contribution from semimuonic decays of D meson pairs, although it 
is significantly smaller. Figure [1] shows the opposite-sign muon pair distribution col- 
lected in p-Pb collisions, compared to the superposition of the different contributions, 
with yields obtained through a multi-step fit procedure, explained below. 

The J /if}, if}', Drell-Yan and open charm contribution shapes have been obtained 
by Monte-Carlo simulation, using a detailed description of the experimental setup 
and of the effects induced by multiple scattering, energy loss, trigger requirements, 
etc. The simulated events were reconstructed with the same programs as the real 
data and had to survive the same selection criteria. 

The Drell-Yan and open charm contributions were generated with the Monte- 
Carlo event generator PYTHIA [TB] (version 6.125) using the GRV 94 LO set of 
parton distribution functions [16] and a primordial k<i Gaussian distribution, with 
width 0.8 GeV/c and 1.0 GeV/c for Drell-Yan and open charm, respectively [17] . 

The J /if) and if>' events were generated with a uniform coshes distribution and 
with Gaussian rapidity distributions, of mean y = —0.2 and width a = 0.85, com- 
mon to all targets. These numerical values were obtained from the analysis of proton- 
nucleus data samples previously collected by NA50, at 450 GeV, with about 10-15 
times higher statistics (the "HI 96/98" data sets, described in detail in Section 6). 
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Figure 1: Measured opposite-sign dimuon mass spectrum for the p-Pb data sample, 
compared with the superposition of expected contributions. 



The pt distributions were generated according to a superposition of a thermal func- 
tion in transverse mass, miK^mT/T) (good for low pr), with a harder component, 
[l + (pt/Po) 2 ] j to describe the high p T tail. The measured pr distributions were 
well described with T values ranging from 266 MeV (for Be) to 288 MeV (for Pb) and 
a common po = 2.349 GeV/c parameter. Figure [2] compares the measured and simu- 
lated y and pi distributions, for dimuons in the mass region 2.9 < M < 3.3 GeV/c 2 , 
where the J/ip contributes more than 95% of all events. 

The dimuon mass distributions of the J/^> and ■*/>' resonances, obtained separately 
for each target, are completely determined by experimental effects (multiple scatter- 
ing and energy loss) . Their mass shapes are extracted from a fit to the reconstructed 
Monte-Carlo events using a pseudo-Gaussian parameterization, i.e., a Gaussian func- 
tion with a mass-dependent width which allows to properly describe the tails of the 
distorted peak. Given the very high statistics of the measured samples in the mass 
window of the resonances, even a small discrepancy between the simulated shapes 
and the measured data would significantly increase the fit % 2 and, thus, affect the 
fitting results. Therefore, the parameters of the function describing the widths and 
tails of the charmonia mass distributions are eventually left free and adjust on the 
observed real data (see details in Ref. [T3]). 

The Monte-Carlo simulation also provides the detection acceptances of each rel- 
evant physics process, as the ratio between reconstructed and generated events, 
in the kinematical domain mentioned above, which amount to 13.5%, 16.2% and 
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Figure 2: Measured and simulated y (left) andpx (right) distributions for J/tp events. 



14.4%, respectively for the J/V>, if >l and Drell-Yan dimuons (in the mass window 
2.9-4.5 GeV/c 2 ). The Monte Carlo generation functions used to compute the char- 
monia acceptances were tuned in such a way that they describe the experimental 
data (see again Fig. [2]). Since our final results have shown to be essentially insensi- 
tive to reasonable small changes of the Monte-Carlo distributions, no other systematic 
uncertainty has been assigned to the acceptance calculation apart from the one aris- 
ing from the finite generated sample, reported in Table [H The Drell-Yan acceptances 
change by ~0.5-1.5% in the 2.9-4.5 GeV/c 2 mass region, systematically for all targets, 
when computed using PYTHIA and different parameterizations of parton distribution 
functions. 

The opposite-sign dimuon mass spectra were fitted, independently for each data 
sample, to the superposition of the contributions previously mentioned, in the mass 
region between 1.5 and 8.0 GeV/c 2 , using the maximum likelihood method. The com- 
binatorial background and the empty target components were kept fixed, while the 
normalisations of the "signal contributions" (J/V, Drell-Yan and open charm) were 
left as free parameters. The systematic uncertainties in the extraction of the number 
of events were investigated fitting the mass spectra under different approaches, such 
as various descriptions of the mass range below 2.3 GeV/c 2 , changing the final fit 
mass range, checking the influence of reasonable parameterizations of J/ip tails under 
the continuum, using different parton distribution functions for Drell-Yan genera- 
tion, etc. From the previous study, we were able to select the best inputs for the 
fit procedure based on the stability of the results as well as on a x 2 / n df criterium. 
However, the xV n df insensitiveness to the used parton distribution functions and to 
reasonable different tails of the J/ip resonance led to somewhat arbitrary choices. As 
a consequence, systematic uncertainties due to these sources must be assigned to the 
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measured number of events, as described in Section 4 and numerically presented in 
Table El The final adopted analysis procedure provides the best overall description of 
the experimental data. Table [3] reports the number of J/ip, ip' an d Drell-Yan events 
(with their statistical error) fitted from each data set, together with the corresponding 
% 2 /ndf values showing the good quality of the fits. 
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Table 3: Number of J/ip, ip' and Drell-Yan events, with their statistical error, and 
reduced x 2 °f eacn data sample fit. 



4 Cross-sections results 

The production cross-sections are derived from the fitted number of events, the in- 
tegrated luminosity, the detection acceptances and from the different efficiencies af- 
fecting our measurement and analysis procedure, already discussed and presented in 
the previous sections. 

Table |4] shows the resulting J/ip, if > ' and Drell-Yan production cross-sections, per 
target nucleon, in the phase space window of our measurement. The values for char- 
monia include the branching ratios of their decay into muons. The Drell-Yan produc- 
tion cross-sections presented in this table (given for the mass range 2.9-4.5 GeV/c 2 ) 
have been multiplied by "isospin correction factors" [6] between 0.954 and 0.947. 
The isospin corrected cross-sections are equivalents to the ones that would have been 
measured if the target nuclei were made of protons only. 

Since the data were collected with relatively frequent target changes, most time- 
dependent systematic effects (such as the ones affecting the beam intensity measure- 
ment) should be identical among the six data sets. With respect to the data samples 
previously collected by NA50, when each p-A system was studied in a different data 
taking period (often meaning different years), this running mode significantly reduces 
the systematic uncertainty in the determination of the nuclear dependences, where 
only the relative changes of the production yields, from target to target, are impor- 
tant. This is the reason why the systematic errors quoted in Table H] only refer to 
uncertainties specific to each data set. These uncertainties may be separated in two 
different kinds of sources: the ones originated by the experimental setup performance 
(~ 0.4-0.7% from the target material properties; ~0.2% from the reconstruction ef- 
ficiencies; ~2.3% from the dimuon trigger efficiency), and the ones related to the 
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Table 4: Production cross-sections, per target nucleon, for J/ijj and ip' production 
(times branching ratios into muons), and for Drell-Yan dimuons in the mass range 
2.9-4.5 GeV/c 2 , integrated in the phase space window of our measurement (—0.425 < 
y cm < 0.575 and — 0.5 < coshes < 0.5). The values in parentheses indicate the relative 
systematic uncertainties which are not common to all targets. 

extraction of the physical yields, as described in the Section 3 (0.2% "resp. 0.4%" 
for the charmonia "resp. Drell-Yan" acceptance and ~ 0.4-1.8 % "resp. <1 %" for ip' 
"resp. Drell-Yan" due to the J/ip high mass tail uncertainty). Still concerning the 
Drell-Yan yields, their low statistics do not allow us to choose which set of parton 
distribution functions better describes the experimental results. This gives rise to 
a -0.3% to +2.8% uncertainty (including acceptance corrections), to be considered 
when comparing our results with those obtained from other data sets, and which 
affects systematically changes the results of each and every target in the same way. 

Finally, the luminosity was computed from the integrated number of incident 
protons within the spectrometer lifetime (see Section 2 for details) and from the 
effective number of target particles per unit surface. The uncertainty on the number 
of incoming protons was estimated to be 3 % based on the systematic error of the 
latest calibration of the ionization chambers |13j . This value also affects all targets 
in the same way and was not included in Table HI However, it must be taken into 
account when comparing these results with those obtained from data collected in 
previous years, in what concerns their absolute normalisation. 

All the different sources specific to each of the targets are summarized in Table [5] 
and a detailed discussion can be found in Ref. [T3] . 

Table [6] presents the ratios between the J/ip, ip' and Drell-Yan production cross- 
sections. Their uncertainties are dominated by the statistical errors of the Drell-Yan 
and ip' values. 

5 Nuclear dependences 

The main topic of this work is the study of the nuclear dependence of the J/ip and ip' 
production cross-sections. It is generally accepted that charmonia production, essen- 
tially due to gluon fusion, is a hard process and, therefore, is basically independent of 
the amount of nuclear matter surrounding the production point. This is also the case 
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Table 5: Contributions to the overall systematic uncertainty of cross-section mea- 
surements, specific to each of the targets. 



CT (DY2.9-4.5) q(DY 2 .9-4.5) B^<t(VQ 

p-Be 54.6 ± 3.2 (0.5%) 0.953 ± 0.082 (0.5%) 1.745 ± 0.085 (0.6%) 

p-Al 52.5 ± 2.8 (0.6%) 0.991 ± 0.076 (0.5%) 1.889 ± 0.078 (0.7%) 

p-Cu 49.3 ± 2.5 (0.9%) 0.786 ± 0.063 (1.1%) 1.593 ± 0.076 (1.9%) 

p-Ag 46.0 ± 2.4 (0.5%) 0.735 ± 0.061 (0.6%) 1.599 ± 0.084 (0.7%) 

p-W 44.9 ± 2.1 (0.6%) 0.639 ± 0.052 (1.1%) 1.422 ± 0.076 (1.5%) 

p-Pb 47.0 ± 1.9 (0.7%) 0.687 ± 0.047 (1.2%) 1.461 ± 0.064 (1.6%) 



Table 6: Ratios of J/V, ip' and Drell-Yan production cross-sections, in our phase 
space window (—0.425 < y cm < 0.575 and —0.5 < coshes < 0.5). The values in 
parentheses indicate the relative systematic uncertainties. 

of high mass dimuon (Drell-Yan) production, where a linear scaling of the produc- 
tion cross-section with the target nucleus mass number, in p-A collisions, has been 
observed [21 [18] . This linearity is usually expressed by saying that a is unity, where 
a is the exponent in the power law o v a = &oA a (the "a parameterization"). It has 
been observed that a decreases as a function of xf and increases as a function of 
Pt (the "Cronin effect"). Since these changes affect Drell-Yan production (where no 
final state effects are expected), they are attributed to initial state effects: nuclear 
modifications of the parton distribution functions (shadowing, etc.) and energy loss 
of the incoming partons as they traverse the target nucleus [Sj. Charmonia produc- 
tion, besides these initial state effects, is also sensitive to nuclear absorption of the 
produced cc bound state, before and after forming the final J/ip or if)' resonances, 
through strong interactions with the nuclear environment on their way out of the 
target nucleus. 

The left panel of Fig. [3] shows the Drell-Yan production cross-section, per target 
nucleon, as a function of the target mass number, as extracted from our six data 
samples. The flat pattern confirms the linear scaling previously observed in several 



9 



proton-nucleus experiments [21 EE], including NA50 [B], indicating that initial state 
effects play a negligible role in Drell-Yan production at our energies and in the NA50 
phase space window (mid-rapidity). From the present analysis we extract q?dy = 
0.982 ± 0.021. The right panel of Fig. [3] shows the corresponding results for the 
J /ip and if)' production cross-sections, which lead to a values smaller than unity: 
a J/lP = 0.925 ± 0.009 and ay = 0.852 ± 0.019. 
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Figure 3: Drell-Yan (left) and charmonia (right) production cross-sections, per nu- 
cleon, as a function of the target mass number. Lines are fits to the a parameterization 
or to the Glauber model, as indicated. 



A much more accurate formalism to describe the observed charmonium produc- 
tion cross-sections is provided by the Glauber model [19] . Within this framework, a 
proton-nucleus or nucleus-nucleus interaction is considered as a set of independent 
interactions between the projectile and target nucleons, assuming that the properties 
of the nucleons do not change after the first collision and that they can interact fur- 
ther with the same cross-section. The basic input for the Glauber formalism in p-A 
collisions is the target nuclear thickness, T^(b), representing the number of nucleons 
per unit surface at impact parameter b. The nuclear density profiles of light nuclei 
(A < 17) were parameterized using the Fermi oscillator model. For the other nuclei 
we used Woods-Saxon distributions, with the parameters given in [20], except for 
the Ag nucleus, missing in those tables, for which we have used the parameterization 
of the FRITIOF 7.02 Monte Carlo code ED. The data were then fitted with the 



expression 



a 
~A 



(A-r 



G 
abs 



db(l-eM-(A-l)^ hs T A (b))) 



leaving as free parameters the normalisations, a£, and the interaction cross-sections 



of the or ip' on their way through the nuclear matter, a^ hs 
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The resulting fits to the J/ip and ip' production cross-sections are presented in 
Figs. [3] (right) and HI respectively as a function of A, the target mass number, and 
L, the average distance of nuclear matter crossed by the charmonium states on their 
way out of the nucleus. The solid lines join together the results of the Glauber 
calculations, made independently for each individual target, while the dotted lines 
represent the error bands resulting from the uncertainties on the normalisations and 
absorption cross-sections. 

The detailed Glaubel model formalism can be approximated, in first order, by the 
"pL parameterization" , 



a 



pA 



pL 



A exp(- 



-o 



pL 

abs 



(PL)) 



where a^ bs is the (approximate) absorption cross-section of the charmonia state and 
p the density of the nuclear matter. L is computed, for each target nucleus, in the 
framework of the Glauber model, as explained in [7]. The corresponding fits, shown in 
Fig. H] as dashed lines, are almost indistinguishable from the full Glauber calculations, 
contrary to what happens with the a parameterization, which significantly diverges 
from the correct calculation in the region of very light nuclei, as shown in Fig. [3] (right) 
and discussed in detail in [7J, [33] • The results of the fits are collected in Table [7J 
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Figure 4: Measured J/ip and production cross-sections, per target nucleon, fitted 
with the full Glauber model and with the pL approximation. 



We see that the if)' suffers a stronger nuclear absorption than the J/ip. However, 
the relative differences between the nuclear dependences of the two charmonium states 
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a fit 


erg a X 2 /ndf 


if/ 


5.6 ± 0.2 nb/nuc 0.925 ± 0.009 0.8 
123 ±q0 pb/nuc 0.852 ± 0.019 2.0 




Glauber fit 


^o G ^ (mb) xVndf 


J/V> 
if/ 


5.1 ± 0.1 nb/nuc 4.6 ± 0.6 0.4 
101 ± 6 pb/nuc 10.1 ± 1.6 1.7 




pL fit 


<L (mb) xVndf 


J/V 


5.0 ± 0.1 nb/nuc 4.1 ± 0.5 0.3 
99 ± 5 pb/nuc 8.2±1.0 1.5 



Table 7: Charmonia nuclear dependence parameters, extracted with the a- 
parameterization, the Glauber model, and its pL approximation. 



can be studied with better accuracy through the ratio B'^a(ip')/B^a(J/ip), shown 
in Fig. [5] as a function of A (left) and L (right). This ratio cancels out part of the 
systematic uncertainties, which affect both charmonium states in the same way. 
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Figure 5: Ratio between the ip' and J/ip production cross-sections, as a function of 
A (left) and L (right). 

Fitting these measurements as a function of the target mass number, with the 
a parameterization, or as a function of L, with the pL parameterization, we can 
directly obtain the parameters which characterize the relative nuclear dependences: 
Aa = a(if/) — a(J/if>) and Acrf bg = <J^ s (ip') — crf bs (J/^>). ^ i s n °t possible to directly 
extract Acr^ s with a Glauber fit. The corresponding results are collected in Table El 
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a fit 


cW/JAA) (%) 


Aa 






2.21 ±0.17 


-0.076 ±0.018 


1.9 




pL fit 


(%) 


A^at (mb) 


X 2 /ndf 




1.98 ±0.10 


4.2 ± 1.0 


1.7 



Table 8: Relative absorption measured directly from the B' cr (if)') /B w <r (if) ratios, 
and fitted with the a and pL parameterizations. 

The nuclear absorption of the charmonium states can also be studied through the 
ratios of their production cross-sections with respect to the Drell-Yan cross-section, 
J /if) /DY and if)' /DY. Under the assumption (supported by existing results [2J [6]) 
that Drell-Yan production has negligible initial and final state effects, the values ex- 
tracted from these ratios should be essentially identical to those obtained from the 
J /if) and ip' absolute production cross-sections. The rather low Drell-Yan statistics of 
the data samples collected in year 2000, however, do not allow us to make significant 
comparisons between the two approaches. For instance, the Glauber formalism ap- 
plied to the J/if> / DY ratios gives <J^ hs (J/ip) = 3.4 ± 1.2 mb, while the absolute J /if) 
cross-section analysis gives 4.6 ± 0.6 mb. Nevertheless, for completeness, we give the 
values extracted from the analysis of the J/if> / DY and if)' / DY ratios at the end of 
the next section (Table [T3|) . together with the values extracted from the simultaneous 
analysis of all NA50 p-A data sets. 

6 Comparison with previous NA50 p-A results 

We will now compare the results shown in the previous sections with those obtained 
by NA50 using proton-nucleus data collected in previous years. For the purpose of 
this comparison, we will label the data sets analysed in detail in the present paper 
as "HI 2000" (given the "high intensity" of the proton beam and the year of data 
taking). 

NA50 collected its first p-A data sets between 1996 and 1998, with a high intensity 
450 GeV proton beam. Five different targets were used (Be, Al, Cu, Ag and W), each 
one during a period of around one week. These data samples, referred as "HI 96/98", 
are the ones with the highest statistics (between 10 and 15 times higher than the 
data sets analysed in the present paper). 

The high intensity beam induced high hit multiplicities in the wire chambers, 
which the standard offline reconstruction software used by NA50 in the first years 
had difficulties in handling properly. This motivated a significant upgrade of the 
reconstruction software, after which the data could be reconstructed with higher 
efficiency. The analysis details and results can be found in [7J [T4] . 

In parallel, new p-A data samples were collected at a reduced beam intensity 
(~ 1-5 10 s protons per burst). The energy of the beam was 450 GeV and the targets 
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used were the same as before. These data sets, referred hereafter as "LI 98/00", 
have 20-30 % of the statistics of the previous ones. Two independent analyses were 
performed using these data, in slightly different kinematical windows, and the results 
have been published in [6l [7] . 

Given the improvements in the reconstruction software, the "HI 2000" data sets 
were collected at high beam intensities, even slightly higher than those of the "HI 96/98" 
samples. However, the short beam time available prevented the collection of sizeable 
statistics. On the other hand, the use of all the targets in the same data taking 
period significantly contributes to a decrease in the systematic uncertainties, when 
determining the nuclear dependences of the production cross-sections. 

Given the slightly different analysis procedures used, the three data samples can- 
not be directly compared without first performing a few corrections. There are three 
main differences among the "HI 96/98", "LI 98/00" and "HI 2000" data sets or anal- 
yses. First, the rapidity window considered for the extraction of the 450 GeV results 
is either —0.4 < y cm < 0.6 or —0.5 < y cm < 0.5 (we remind that the 400 GeV probe 
the range —0.425 < y cm < 0.575). Second, the Drell-Yan calculations were done 
with different parton distribution functions: MRS G [22] , MRS A (Low Q 2 ) [23] and 
GRV 94 LO [16]. Third, the Monte Carlo simulations of the charmonium states used 
somewhat different y and pt distributions. 

For a coherent comparison of all the three data sets, we corrected the 450 GeV 
results to a common rapidity domain, —0.5 < y cm < 0.5, and also so that they corre- 
spond to the conditions used in the "HI 2000" analysis, as described in the previous 
sections of this paper: GRV 94 LO parton distribution functions and improved char- 
monia rapidity and pt distributions. These corrections, described in detail in [13] , are 
rather small: 0.9 % (0.2 %) for the charmonia (Drell-Yan) y window, between —2.0 % 
and +1.3 % for the Drell-Yan shape imposed by the different PDFs and ~ 1 % for the 
charmonia kinematical distributions. Finally, the HI 96/98 Drell-Yan results were 
corrected by 6.9±1.6%, in order to account for their systematically lower normalisa- 
tion with respect to the other data sets, as reported in [T4"] . 

The change in the function used to generate the rapidity distribution of the 
charmonium states is worth being explained in some detail. The analyses of the 
"HI 96/98" and "LI 98/00" data sets were done assuming that the J/ip mesons were 
produced with a centre of mass rapidity distributed according to a Gaussian centred 
at mid-rapidity and with a r.m.s. of 0.75 [H] or 0.60 [23], independently of the target. 
After these analyses were published [61 [7], we tried to understand better the specific 
issue of the rapidity distributions, from a more detailed look at our measurements, 
despite the fact that our rapidity coverage is relatively narrow (only one unit in ra- 
pidity). Figure [6] shows the rapidity distributions (corrected for acceptance) of events 
in the mass range 2.86 < M < 3.34 GeV/c 2 , from the "HI 96/98" (highest statistics) 
data sample. These distributions were fitted to Gaussians with free widths and either 
centred at mid-rapidity (yo = 0, dashed lines) or with a free mean (solid lines). All 
the five different targets, from Be to W, systematically prefer a significantly nega- 
tive mean value, y ~ —0.2, and a common width, o = 0.85, with reasonable fit 
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quality, x 2 / n df ~ 1-3. On the contrary, forcing the y distributions to be centered 
at mid-rapidity leads to a bad description of the measured data, with x 2 /ndf values 
between 18 and 51, depending on the target. The Monte Carlo simulation done in 
the present work (for the J/ip and the if)') used Gaussian rapidity distributions with 
mean —0.2 and a = 0.85. Nevertheless, we have seen that the final results are es- 
sentially insensitive to reasonable changes in the rapidity distributions used for the 
event generation. 
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Figure 6: J /if; rapidity distributions, in the centre of mass frame, as measured by 
NA50 through the "HI 96/98" data sets, with five different nuclear targets and a 
450 GeV proton beam. The data points are fitted to Gaussian functions with a free 
width and a mean which is either left free (solid lines) or fixed to (dashed lines). 

Table [9] compiles the "HI 96/98" and "LI 98/00" results, after the corrections just 
mentioned. The "LI 98/00" results shown here are the average of the two indepen- 
dent analyses mentioned above, since they are in agreement with each other, when 
compared in a coherent way. 

Table [10] shows the results of each of the three data sets, independently, in what 
concerns the nuclear dependences of the J /if) and if)' production cross-sections, as 
determined using the detailed Glauber formalism. It also shows, in the last line, 
the results obtained from a simultaneous fit to the three data sets, with a common 
absorption cross-section, a^ hs , and two different normalisations, accounting for the 
different energies and rapidity domains (we remind that the 400 GeV data is anal- 
ysed in the window —0.425 < y cm < 0.575). The "HI 2000" absolute cross-section 
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(nb/nucleon) 


(pb/nucleon) 


a(DY 2 . 9 _ 4 . 5 )/A 
(pb/nucleon) 




HT 96/98 T,T 98/00 

111 c/U/ c/O iJl C/O/ UU 


HT 96/98 T,T Q8/DD 

111 c/U/ C/O iJl C/O/ WW 


HT 96/98 T,T 98/flfl 

111 C/U/ C/O iJl C/O/ uu 


Be 
Al 

Cu 
Ag 
W 


5.11±0.18 5.27±0.23 
4.88±0.23 5.14±0.21 
4.74±0.18 4.97±0.22 
4.45±0.15 4.52±0.20 
4.05±0.15 4.17±0.37 


88.2±3.8 89.4±5.3 
84.3±4.6 88.0±4.8 
77.7±3.4 82.1±4.9 
69.8±2.7 72.8±4.3 
61.7±2.6 66.3±7.9 


93.3±4.1 96.6±5.6 
96.1±5.0 104.0±5.9 
101.3±4.6 97.3±6.8 
94.3±3.7 100.4±5.3 
92.7±4.0 96.3±8.9 



Table 9: Updated J/^, ip' an d Drell-Yan production cross-sections, in p-nucleus 
collisions at 450 GeV, as extracted from the "HI 96/98" and "LI 98/00" data sets, 
in the —0.5 < y cm < 0.5 and —0.5 < coshes < 0.5 phase space window. 



systematic errors, affecting all targets in the same way, are not considered in the 
global fit, since the 400 GeV data has its own free normalisation. 





a G (J/V) 
(nb/nuc) 


^bs(J/V>) 

(mb) 


X 2 /ndf 


(pb/nuc) 


(mb) 


X 2 /ndf 


HI 96/98 


5.5 ± 0.2 


4.4 ± 1.0 


0.9 


101 ± 5 


7.6 ± 1.4 


1.3 


LI 98/00 


5.7 ± 0.3 


4.1 ± 1.4 


0.6 


100 ± 7 


5.7 ± 2.0 


0.6 


HI 2000 


5.1 ± 0.1 


4.6 ± 0.6 


0.4 


101 ± 6 


10.1 ± 1.6 


1.7 


450 GeV 
400 GeV 


5.6 ± 0.1 
5.1 ± 0.1 


4.5 ± 0.5 


0.7 


106 ± 4 
95 ± 4 


8.3 ± 0.9 


1.3 



Table 10: Nuclear dependence parameters of the J/ip and ip' production cross- 
sections, from each of the three data sets and from a global fit to all of them. 

Figured shows together all the J/ijj (left) and ip' (right) absolute production cross- 
sections measured by NA50 in p-nucleus collisions, as collected in Tables H] and EE The 
curves are the result of the simultaneous Glauber fits to all the data points and the 
error bands represent the uncertainties in the absorption cross-section and in the 
normalisations (see Table [TU1 for the numerical values). 

If instead of using the full Glauber calculation we use the approximate pL param- 
eterization, a does not change but the cr^ bs values are smaller, by ~ 10 % for the J ftp 
and by ~ 15 % for the ip'. 

The pL parameterization is particularly useful to determine the relative absorption 
cross-section of the ip' with respect to that of the J/ip, Acx^,. In the same way, the a 
parameterization is also used to quantify the relative charmonia absorption through 
the study of the Aa parameter. 

Table dU shows the Acr£ bs and Aa values independently obtained from each of 
the three data sets and also when we perform a simultaneous fit to all the data, 
assuming that the ratio between the ip' and the J/ip production cross-sections does 
not change from 400 to 450 GeV. We can see that the values given by the global fit, 
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Figure 7: J ftp (left) and ip' (right) cross-sections from the three different NA50 data 
sets. The full line represents the result of a simultaneous Glauber fit with a common 
absorption cross-section and two independent normalisations. The error bands reflect 
the absorption cross-section and normalisation uncertainties. 









X 2 /ndf 


^ WW 


Aa 


X 2 /ndf 




(%) 


(mb) 




(%) 






HI 96/98 


1.83±0.06 


2.6±0.6 


0.5 


1.95±0.09 


-0.045±0.011 


0.7 


LI 98/00 


1.79±0.08 


1.9±1.0 


0.4 


1.87±0.13 


-0.032±0.017 


0.5 


HI 2000 


1.98±0.10 


4.2±1.0 


1.7 


2.21±0.17 


-0.076±0.018 


1.9 


Global 


1.85±0.04 


2.8±0.5 


0.9 


1.98±0.07 


-0.048±0.008 


1.0 



Table 11: Relative nuclear absorption of the tp' with respect to the J/ip, obtained 
from the B' a(ip') /B^afy) ratios, using the pL and a parameterizations. 



Acre's = 2.8 ± 0.5 mb and Aa = —0.048 ± 0.008, are essentially the same as the ones 
determined from the high statistics "HI 96/98" data set, on its own. 

From the global fits leading to Table EH we see that the ratio between the char- 
monia production cross-sections at 400 GeV, in the —0.425 < y cm < 0.575 window, 
and at 450 GeV, in the —0.5 < y cm < 0.5 window, is 0.90 (with a 3% relative error 
bar). This value is in perfect agreement with the expectation of a NLO QCD cal- 
culation made in the Color Evaporation Model framework [25]. Using the "Schuler 
parameterization" [26], which provides a phenomenological description of J/ip pro- 
duction cross-sections as a function of yfs, we would expect 0.933 ± 0.008, a value 
which includes the corresponding changes in the x-p distributions and that is quite 
close to the result derived from our analysis. 
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As already briefly mentioned at the end of the previous section, it is interesting to 
study the J/ip and ip' nuclear dependences using the ratios of production cross-sections 
B^a(J/ip)/a(DY 2 . 9 ^. 5 ) and B'^a(ip')/a(DY 2 . 9 _ 4 . 5 ). In what concerns the "HI 2000" 
data set, these values have already been presented in Table [6j The corresponding 
values for the two other data sets are collected in Table [T2j They were obtained from 
the analyses described in the references [21 El EH] and corrected, where appropriate, 
as explained above. 





-Bum " 


(J/tf) 






*m 




a(DY 2 


9-4.5) 


a(DY 2 


9-4.5) 




'(JAM 




HI 96/98 


LI 98/00 


HI 96/98 


LI 98/00 


HI 96/98 


LI 98/00 


Be 


54.7±1.5 


54.6±1.8 


0.946±0.039 


0.932±0.051 


1.73±0.05 


1.70±0.07 


Al 


50.8±1.2 


49.9±1.7 


0.876±0.030 


0.859±0.044 


1.73±0.04 


1.71±0.07 


Cu 


46.7±1.0 


51.0±2.0 


0.768±0.024 


0.848±0.044 


1.64±0.03 


1.66±0.05 


Ag 


47.2±1.0 


45.5±1.3 


0.740±0.021 


0.735±0.034 


1.57±0.03 


1.61±0.05 


W 


43.7±0.9 


43.6±1.7 


0.666±0.020 


0.665±0.044 


1.53±0.04 


1.52±0.07 



Table 12: Ratios of J/tp, ip' and Drell-Yan production cross-sections, in p-A collisions 
at 450 GeV, as extracted from the "HI 96/98" and "LI 98/00" NA50 data sets, in 
the —0.5 < y cm < 0.5 and —0.5 < cos 9 C s < 0.5 kinematical window. 

From the cross-section ratios of Tables [6] and d2] we can extract the nuclear ab- 
sorption parameters, using the Glauber formalism, as we have done from the absolute 
cross-sections (in Table [TO]) . The results are summarised in Table [TBI where the last 
line shows the values resulting from a global fit to the three data sets. The derivation 
of the charmonia nuclear absorption parameters using the ratios of cross-sections, 
Bfj,^cr(J /ip) I '<7(-D>2.9-4.5) and B'<t(i(/) j '<r(-DY"2. 9-4.5), is very important to establish 
the "normal nuclear absorption curve" in heavy-ion collisions, where the suppression 
pattern of the J/ip and if)' resonances is determined, for each centrality bin, using the 
ratio of their production yields to the corresponding yield of Drell-Yan events. 





^(J/V/DY) 




xVridf 


a G (^'/DY) 




X 2 /ridf 




(nb/nuc) 


(mb) 




(nb/nuc) 


(mb) 




HI 96/98 


57.7 ± 1.7 


4.3 ± 0.7 


1.5 


1.07 ± 0.05 


7.8 ± 1.2 


0.8 


LI 98/00 


58.2 ± 2.3 


4.4 ± 1.0 


0.9 


1.05 ± 0.07 


6.8 ± 1.7 


0.8 


HI 2000 


57.2 ± 3.4 


3.4 ± 1.2 


0.4 


1.14 ± 0.11 


8.8 ± 2.3 


0.9 


450 GeV 
400 GeV 


57.5 ± 1.3 
59.3 ± 1.8 


4.2 ± 0.5 


0.7 


1.07 ± 0.04 
1.09 ± 0.05 


7.7 ± 0.9 


0.7 



Table 13: Nuclear dependence parameters of the J/ip and ip' resonances, using cross- 
section ratios, from each of the three data sets and from a global fit to all of them. 

It is important to note that the J/ip absorption cross-section derived from the 
J/ijj /DY ratios, c^ )S (J/'0) = 4.2 ± 0.5 mb, is in excellent agreement with the value 
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derived from the absolute production cross-sections, 4.5 ± 0.5 mb (Table [TUT) . Also 
the if)' values obtained in the two alternative approaches are in very good agreement, 
ag s (ip') = 7.7 ± 0.9 mb (Table [13]) and 8.3 ± 0.9 mb (Table HDP . When defining the 
normal nuclear absorption curve, as a function of centrality, in heavy-ion collisions, 
to compare with the J /if) suppression pattern obtained in those collisions through the 
J /if) /DY cross-section ratio, the value 4.2 ± 0.5 mb should be used, for consistency. 
This a^ hs value is in agreement with previous published results [5J [27J but it should 
be noticed that when NA51 pp and pD charmonia measurements [28] were included 
in the Glauber fit [27], a 10% systematic error had to be added due to the different 
analysis techniques used in NA51 data. The inclusion of the pp and pD values 
in the Glauber fit has a strong influence on the derived value of a^ hs and since a 
coherent analysis including NA51 points cannot be performed any more due to raw 
data unavailability, these results were excluded from the analysis reported here. 

As a final remark, the a^ bs energy and kinematical dependence was also inves- 
tigated using other p-A production cross-sections, published by experiments with 
similar spectrometers: NA3 p[] and NA38 [3J, at = 200 GeV. These values, fit- 
ted with the full Glauber model, lead to cr^ hs (J/ip) = 3.3 ± 3.0 mb, which, within 
the present precision, is compatible with our present measurement. The inclusion of 
these points in the global Glauber fit using J /if) cross-sections does not change the 
final a^ hs value. Nevertheless, the energy dependence of the absorption cross-section 
remains, today, an open question. 

7 Summary and conclusions 

This paper presents the J /if), if)' and Drell-Yan production cross-sections in proton- 
nucleus collisions, as extracted from data collected by NA50 in year 2000, with 
400 GeV protons and six different targets: Be, Al, Cu, Ag, W and Pb. The fre- 
quent rotation of the target material facing the beam significantly decreases the 
time-dependent systematic uncertainties related to beam stability, trigger system, 
etc. The accumulated statistics, however, suffers from the very short beam time 
available, which affects in particular the if)' and Drell-Yan production results. 

The nuclear dependence of the charmonium states production is studied, as a 
function of the target nucleus, with three different descriptions: the a parameteriza- 
tion, the detailed Glauber formalism and its "pL" approximation. The J /if) and if>' 
absorption cross-sections obtained with the Glauber model, from the 400 GeV data 
sets, are: of bs (J/V) = 4.6 ± 0.6 mb and cr® bs (ip') = 10.0 ± 1.5 mb. 

A detailed comparison with previous NA50 p-A data samples, collected at 450 GeV 
with five different nuclear targets, shows that all data can be described with com- 
mon absorption cross-sections. The "HI 2000" data set is particularly important 
in the determination of cr a bs(J/V ; ) from the study of the absolute production cross- 
sections, given the lower systematic uncertainties on the beam normalisations. On 
the other hand, the two data sets previously collected by NA50 have more statistics 
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and, therefore, lead to more accurate measurements of the B^a(J /ip) / a(DY 2 ,g-4^) 
and B' cr(t/j') /<j(DY 2 . 9-4.5) cross-section ratios, which are insensitive to the beam nor- 
malisation. Global fits to the cross-section values give: 0"!^ s (J/-?/>) = 4.5 ±0.5 mb and 
c^ bs (ip') = 8.3 ±0.9 mb. The corresponding results from the B^cr(J /ip) / a(_DY"2.9-4.5) 
and /o"(£ ) 1 / 2.9-4.5) cross-section ratios are in excellent agreement: o~^ hs (J/ip) = 

4.2 ± 0.5 mb and cr^ hs (ip') = 7.7 ± 0.9 mb. The fact that the two alternative ways of 
deriving the charmonia absorption cross-sections, through absolute production cross- 
sections or through their ratios, are in such good agreement (furthermore, being 
mostly determined from different data sets), is an independent and accurate indi- 
cation that, at our energies and in our kinematical phase space window, Drell-Yan 
production is not sensitive to nuclear effects, including initial state effects such as nu- 
clear modifications of the quark distribution functions and energy loss of the incident 
nucleons. It would be interesting to see if a similar study made as a function of the 
dimuon pr or xp would modify this observation. 

For consistency reasons, the values o~^ s (J/ip) = 4.2 ± 0.5 mb and cr^ hs (ip') — 
7.7 ±0.9 mb should be used in the comparison with the cross-section ratios measured 
in heavy ion collisions. A significantly stronger absorption of the ip' resonance with 
respect to the J/ip is seen through the study of the relative production yields between 
the two charmonium states, in our phase space window. 

The results presented in this paper were obtained assuming that final state nuclear 
absorption of the charmonium states, as calculated through detailed Glauber calcu- 
lations with a single free parameter, a^ hs , is the only nuclear effect influencing the 
production of these resonances in p-nucleus collisions. If future measurements (for in- 
stance, of open charm production cross-sections in p-nucleus collisions) demonstrate 
the existence of nuclear effects (such as anti-shadowing) on the gluon distribution 
functions at SPS energies, the numerical values given in this paper should be seen 
as effective values, resulting from the convolution of initial and final state nuclear 
effects. 
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